Search for: All records

Currently, it is challenging to investigate aneurismal hemodynamics based on current in vivo data such as Magnetic Resonance Imaging or Computed Tomography due to the limitations in both spatial and temporal resolutions. In this work, we investigate the use of modal analysis at various resolutions to examine its usefulness for analyzing blood flows in brain aneurysms. Two variants of Dynamic Mode Decomposition (DMD): (i) Hankel-DMD; and (ii) Optimized-DMD, are used to extract the time-dependent dynamics of blood flows during one cardiac cycle. First, high-resolution hemodynamic data in patient-specific aneurysms are obtained using Computational Fluid Dynamics. Second, the dynamics modes, along with their spatial amplitudes and temporal magnitudes are calculated using the DMD analysis. Third, an examination of DMD analyses using a range of spatial and temporal resolutions of hemodynamic data to validate the applicability of DMD for low-resolution data, similar to ones in clinical practices. Our results show that DMD is able to characterize the inflow jet dynamics by separating large-scale structures and flow instabilities even at low spatial and temporal resolutions. Its robustness in quantifying the flow dynamics using the energy spectrum is demonstrated across different resolutions in all aneurysms in our study population. Our work indicates that DMD can be used for analyzing blood flow patterns of brain aneurysms and is a promising tool to be explored in in vivo. 

Abstract Distribution of bed shear stress is the critical factor in regulating the meandering of single-thread rivers. However, the impact of ice cover on bed shear stress is largely unknown. In this study, we develop a theoretical model of cross-stream momentum balance to examine the distribution of bed shear stresses in ice-covered meandering rivers. To validate the theoretical model, field surveys were carried out in a river reach of the Red River in Fargo, North Dakota. Data monitoring was completed using an Acoustic Doppler Current Profiler to obtain time-averaged velocity profiles. Our theoretical model indicates that an ice covering develops high-shear zones near both the inner and outer banks, which might exacerbate sediment transport and enhance bank erosion. Velocity measurements confirm the results of the proposed model and demonstrate a clear impact of meandering river banks on velocity profiles and secondary flow patterns under ice cover. Based on our results, we hypothesize that ice cover increases turbulent stresses near banks, which in turn lead to the enhancement of the bed shear stress. Our work provides new insights into the impact of ice cover on bed shear stress distribution, which could play an important role in driving sediment-transport processes and the long-term morphodynamic evolution of meandering rivers seasonally covered by ice. 

Abstract We performed the transport of a breast cancer cell (MB231-TGFb) in a microvessel using high-resolution simulations. Using open-source imaging software Slicer3D and Meshmixer, the 3D surface mesh forming the cell membrane was reconstructed from confocal microscopic images. The Dissipative Particle Dynamics method is used to model the cell membrane. The extracellular fluid flow is modeled with the Immersed Boundary Method to solve the governing equations of the blood plasma. The unsteady flow is applied at the inlet of the microchannel with an oscillatory pattern. Our results showed that the extracellular flow patterns are highly dependent on the waveform profile. The oscillatory flow showed the creation of vortices that influence the cellular deformations in the microchannel. These results could have implications on the destination of the cancer cells during transport in physiological flows. 

Given the complexity of human left heart anatomy and valvular structures, the fluid–structure interaction (FSI) simulation of native and prosthetic valves poses a significant challenge for numerical methods. In this review, recent numerical advancements for both fluid and structural solvers for heart valves in patient-specific left hearts are systematically considered, emphasizing the numerical treatments of blood flow and valve surfaces, which are the most critical aspects for accurate simulations. Numerical methods for hemodynamics are considered under both the continuum and discrete (particle) approaches. The numerical treatments for the structural dynamics of aortic/mitral valves and FSI coupling methods between the solid Ωs and fluid domain Ωf are also reviewed. Future work toward more advanced patient-specific simulations is also discussed, including the fusion of high-fidelity simulation within vivo measurements and physics-based digital twining based on data analytics and machine learning techniques. 

Transport of cells in fluid flow plays a critical role in many physiological processes of the human body. Recent developments of in vitro techniques have enabled the understanding of cellular dynamics in laboratory conditions. However, it is challenging to obtain precise characteristics of cellular dynamics using experimental method alone, especially under in vivo conditions. This challenge motivates new developments of computational methods to provide complementary data that experimental techniques are not able to provide. Since there exists a large disparity in spatial and temporal scales in this problem, which requires a large number of cells to be simulated, it is highly desirable to develop an efficient numerical method for the interaction of cells and fluid flows. In this work, a new Fluid-Structure Interaction formulation is proposed based on the use of hybrid continuum-particle approach, which can resolve local dynamics of cells while providing large-scale flow patterns in the vascular vessel. Here, the Dissipative Particle Dynamics (DPD) model for the cellular membrane is used in conjunction with the Immersed Boundary Method (IBM) for the fluid plasma. Our results show that the new formulation is highly efficient in computing the deformation of cells within fluid flow while satisfying the incompressibility constraints of the fluid. We demonstrate that it is possible to couple the DPD with the IBM to simulate the complex dynamics of Red Blood Cells (RBC) such as parachuting. Our key observation is that the proposed coupling enables the simulation of RBC dynamics in realistic arterioles while ensuring the incompressibility constraint for fluid plasma. Therefore, the proposed method allows an accurate estimation of fluid shear stresses on the surface of simulated RBC. Our results suggest that this hybrid methodology can be extended for a variety of cells in physiological conditions.